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Abstract. Two setups with interlocked magnetic flux tubes are used to study the evolution of 
magnetic energy and helicity on magnetohydrodynamical (MHD) systems like plasmas. In one 
setup the initial helicity is zero while in the other it is finite. To see if it is the actual linking or 
merely the helicity content that influences the dynamics of the system we also consider a setup 
with unlinked field lines as well as a field configuration in the shape of a trefoil knot. For helical 
systems the decay of magnetic energy is slowed down by the helicity which decays slowly. It 
turns out that it is the helicity content, rather than the actual linking, that is significant for the 
dynamics. 
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Magnetic helicity has been shown to play an important role in the dynamo process 
([Brandenburg fc Subramanianl 120051) . For periodic systems where helicity is conserved 
simulations have shown that with increasing magnetic Reynolds number ReM the sat- 

1/2 (iBrandenburg fc Doblerl l200ll) . 



uration magnetic field strength decreases like Re M 
This is problematic for astrophysical bodies since for the Sun ReM = 10 9 and galax- 
ies ReM = 10 14 . In order to alleviate this quenching t he magnetic helicity of the small 
scale fields needs to be shed ( Brandenburg et al. 20091) . 

In the active reg ions of the Sun twist ed magnetic field lines have been observed 
(jPevstov et al.lll995l ). Later it was shown ( Leka et al.lll996j) that the magnetic fie ld in 



sunspots gets twisted before it emerges out of the surface. (IManoharan et al.lll996l ) and 
(jCanfield et al.lll999l ) demonstrated that helical structures are more likely to erupt into 
coronal mass ejections. This suggests that the Sun sheds helicity. 

The ma gnetic helicity is related to the mutual linking for two non-intersecting flux 
tubes via |Moffatdll969l ) 



H 



A ■ B dV = 2n<i 



where H is the magnetic helicity, B = V x A is the magnetic field in terms of the vector 
potential A, </>i and 02 are the magnetic fluxes through the tubes and n is the linking 
number. The flux tubes may not have internal twist. In the limit of large ReM H is a 
conserved quantity as well as the linking number. 

In presence of mag netic helicity the magnetic energy decay is constrained via the 
realizability condition ( Moffatt 19691 ) which gives a lower bound for the spectral magnetic 
energy 

M{k) ^ k\H(k)\/2fj. with J M{k) dfc = (B 2 )/2fi , J H{k) dk= (A- B), 

the magnetic permeability (j,q, where (.) denotes volume integra ls. 

In this work we extend earlier work ( Del Sordo et al. 2010l ) where the dynamics of 
interlocked flux rings, with and without helicity, was studied as well as a non-interlocked 
configuration. Here we also study a self-interlocked flux tube in the form of a trefoil knot. 
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Figure 1. The three triple ring configurations 
for the initial time. From left to right: inter- 
locked rings with no helicity, interlocked rings 
with finite helicity and non-interlocked rings 
without helicity. The arrows indicate the di- 
rection of the magne tic field. Adapted from 
ijDel Sordo et al.ll2010t ). 




Figure 2. The initial magnetic field 
configuration for the trefoil knot. 



The three-rings setups consist of three magnetic flux tubes. In two configurations they 
are interlocked where in one the helicity is zero and in the other one it has a finite value, 
as shown in Fig. [I] In the third setup we instead consider unlocked rings. Since the rings 
do not have internal twist, the helicity of this last configuration is zero. We also study 
the evolution of a self- inter locked flux tube having the form of a trefoil knot, with finite 
helicity (Fig. [5]). In this case we have H — 3(f> 2 , so the linking number is n — 3/2. All of 
these setups evolve according to the full resistive equations of MHD for an isothermal 
compressible medium. The Alfven time is used as time unit. 

As a consequence of the rcalizability condition the magnetic energy cannot decay faster 
than the helicity. To check the importance of this restriction we plot the ratio 

H(k) = 2fMoM(k)/k\H(k)\ (0.1) 

in Fig. [3] for k — 4, which is the typical scale of the system. The magnetic energy at 
that scale reaches its allowed minimum value already at t = 10, which confirms the 
importance of the realizability condition in this decay process. The setups with finite H 
show a slower decay than the setups with no helicity (Fig. [3]). The decay of the trefoil 
knot follows approximately the same decay law as the other configuration consisting of 
three rings with finite H . Within the simulation time H decays only to about one half of 




T 



Figure 3. Left: Evolution of the normalized magnetic energy for the trefoil knot (solid/red 
line) compared with various three-ring configurations with n — 2 (dash-dotted line), n = 
(dashed/blue line), and the non-interlocked case (dotted/blue line). Right: H(k) with k = 4 for 
the helical triple ring configuration. 
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Figure 4. Magnetic field lines at 5 Alfven times for the trefoil knot. The colors represent the 
magnitude of the magnetic field. Note that internal twist generation is weak. 

the initial value conserving then the topology. During later times field lines reconnect and 
the hclicity seems to go into internal twist, which is topologically equivalent to linking; 
see Fig. 01 

The slow decay of H conserves the topology of the system. The linking is then even- 
tually transformed into internal twisting during magnetic reconnection. Since both non- 
helical setups evolve similarly we conclude that it is mainly the magnetic helicity and 
not the actual linking which influences the dynamics. The helical trefoil knot evolves in a 
similar manner. This confirms the hypothesis that the decay of interlinked flux structures 
is g overned by magnet ic helicity and that higher-order invariants, advocated for example 



by ( Yeates et al. 201dl ). may not be essential for describing this process. 



In conclusion, we can say that magnetic helicity is decisive in controlling the decay of 
interlocked magnetic flux structures. If the magnetic helicity is zero, resistive decay will 
be fast while with finite magnetic helicity the decay will be slow and the speed of decay 
of magnetic energy depends on the speed at which magnetic helicity decays. This is likely 
an important aspect also in magnetic reconnection problems that has not yet received 
sufficient attention. 



References 

Brandenburg, A. & Subramanian, K. 2005, Phys. Rep. 417, 1 
Brandenburg, A. & Dobler, W. 2001, Astron. Astrophys. 369, 329 
Brandenburg, A., Candelaresi, S. & Chatterjee P. 2009, MNRAS 398, 1414 
Canfield, R. C, Hudson, H. S. & McKenzie, D. E. 1999, Geophys. Res. Lett. 26, 627 
Moffatt H. K. 1969, J. Fluid Mech. 35, 117 

Del Sordo, F., Candelaresi, S. & Brandenburg, A. 2010, Phys. Rev. E 81, 036401 
Gibson, S. E., Fletcher, L., Del Zanna G. & et al. 2002 ApJ 574, 1021 

Leka, K. D., Canfield, R. C, McClymont A. N. & van Driel-Gesztelyi, L. 1996, ApJ 462, 547 
Manoharan, P. K., van Driel-Gesztelyi, L., Pick, M. & Demoulin P. 1996, ApJ 468, L73 
Pevtsov, A. A., Canfield, R. C. & Metcalf, T. R. 1995, ApJ 440, L109 
Yeates, A. R., Hornig, G., & Wilmot-Smith, A. L. 2010, Phys. Rev. Lett., 105, 085002 



